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ABSTRACT 

Coronal mass ejections (CMEs) with angular width > 30° have been ob¬ 
served to occur at a higher rate during solar cycle 24 compared to cycle 23, 
per sunspot number. This result is supported by data from three independent 
databases constructed using Large Angle and Spectrometric Coronagraph Exper¬ 
iment (LASCO) coronagraph images, two employing automated detection tech¬ 
niques and one compiled manually by human observers. According to the two 
databases that cover a larger held of view, the enhanced CME rate actually began 
shortly after the cycle 23 polar held reversal, in 2004, when the polar helds re¬ 
turned with a 40% reduction in strength and interplanetary radial magnetic held 
became ~ 30% weaker. This result is consistent with the link between anoma¬ 
lous CME expansion and heliospheric total pressure decrease recently reported 
by Copalswamy et al. 

Subject headings: Sun: activity. Sun: coronal mass ejections (CMEs), Sun: mag¬ 
netic helds. Sun: corona. Sun: photosphere. Sun: heliosphere 


1. Introduction 


Coronal mass ejections (CMEs) are large, twisted magnetic structures that are ex¬ 
pelled from the solar surface out into the heliosphere at speeds of hundreds of km s“^. 
They are believed to remove magnetic helicity and free magnetic energy from the corona 
that would otherwise b e difficult to eliminate, owing to the high electrical conductivity of 
the corona flLowl 20011) . They cause the most dangerous space weather e hects on Eart h 
( Webb fc Howard 2012 ). First observed in coronagraph images in the 197 0s (Tousev 19^. 
they have since been observed systematically by numerous observatories fjWebb fc Howard 
19941) . incl uding since 1996 the L arge Angle and Spectrometric Coronagraph Experiment 
(LASCO) ( Brueckner et al. 1995 ) on NASA’s Solar and Heliospheric Observatory (SoHO) 
satellite. 
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The relationship between the behavior of coronal mass ejec tions fCMEs) ari d the g lobal 
solar magnetic held has long been a topic of lively interest. ITaihmann et ahl (119981) and 
Li et ahl (120011 ) compared the locations of front-side CMEs observed by LASCO with the 
distribntions of newly open magnetic held lines in the corona, inferring the opening held 
lines from potential-held source-snrface models based on magnetogram observations taken 
before and after each CME. In some cases they fonnd an intrigning morphological similarity 
between the CME emission stru ctnre and the newly-open held lines. Using LASCO data, 
Gopalswamv et ahl (j2003bl 120121) have reported an eqnatorward latitnde ohset between CME 
trajectories and the locations of their associated prominence ernptions, i mplying inhnence on 


the ern ptions from the expanding polar coronal hole helds. Fnrthermore, iGopalswamv et al. 


(j2003aj) fonnd a close relationship at both poles between the cycle 23 magnetic polarity rever¬ 
sals and the statistics of high-latitnde GMEs: the high-latitnde ernptions stopped in Novem¬ 
ber 2000 and May 2002 in the northern and sonthern hemispheres, respectively, ronghly 
coinciding with the north and sonth polar held reversals. 

Several years ago it became clear that the cycle 23/24 minimnm polar helds were approx¬ 
imately 40% weaker than the cycle 22/23 minirnum polar helds, acco mpanied by a 40-45% 

de Tomall201ll). and a 30% d e crease 


decrease in polar coronal hole area ( Wang et ah 20091: 


in the radial interplanetary mean held (IMF ISmith fc Balogh 


20081 


Lnhmann et ahl (120111) 


snggested that the observed enhanced GME rate of cycle 24 may be explained in terms 
of the weak polar helds allowing weaker acti ve region and prominence helds to ernpt, and 
more ejections to escape into the heliosphere. iPetriel (120131) fonnd evidence in the Gompnter 
Aided GME Tracking (GAGTns) and Solar Ernptive Events Data System (SEEDS), GME 
databases that the enhanced GME rates indeed correlated well in time with the weakening of 
the polar helds. In earlier data collected dnring cycle 21 and the rise of cycle 22 (1975-1989) , 
the GME rate was very well correlated with the snnspot nnmber (jWebb fc Howardlll994l) . 
The LASGO GME rates based on the GAGTns and SEEDS databases, in contrast, were sys¬ 
tematically higher per snnspot nnmber after the cycle 23 polar reversal compa red to before , 
and the polar helds decrea sed in strength by abont 40% dnring that reversal (jPetrid 120131) . 
Wang fc Golaninnol (120141) pointed ont that the change in LASGO image cadence beginning 
in 2010, when the image rate approximately donbled, may have cansed the increased rate 
of GME detections. Assnming that the GME detection rate is proportional to the LASGO 
image rate, they normalized the GME detection rate and fonnd a mnch higher correlation 
with the snnspot nnmber, and conclnded that the polar helds have no signihcant ahect on 
the occnrrence rate of GMEs. 


Since then, some evidence has appeared that the cycle 24 GME rate is genninely en¬ 
hanced over the cycle 23 rate, and a physical explanation has been ohered. According to the 
Goordinated Data Analysis Workshop (GDAW) database, the halo GME rate has been higher 
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during cycle 24 than during cycle 23 flGopalswamv et al.N2015l) . even though the sunspot 
number has been smaller by around 40%. The distribution of CME sources in apparent 
longitude has also been much flatter, with proportionally twice as many halo CMEs origi¬ 
nating from central meridian distances > 60°. Howeve r, the average CME spee d and flare 
size have been unchanged compared to cycle 23, leading ICopalswamv et ahl (120151 1 to suggest 
an explanation based on the ambient medium and not source magneti c helds: a decrease 


i n tot al (magnetic -|- plasma) pressure in the corona and heliosphere. iGopalswamv et ah 


(1201411 found evidence for this in a study of CME widths and velocities: the linear correla¬ 
tion between these two quantities has changed little from cycle 23 (Pearson linear correlation 
coefficient r = 0.67 to cycle 24 (r = 0.72), but the constant of proportionality has changed 
by 50%. Cycle 24 CMEs are signihcantly wider than their cycle 23 counterparts without 
being signihcantly faster. Based on NASA OMNI solar wind data taken in situ at 1 AU, 
and assuming that the magnetic held strength, proton density and proton temperature de¬ 
cay with radial distance R from the Sun as R~^, R~^ and respectively, the authors 

argued that the heliospheric total pressure has declined by around 40% between cycles 23 
and 24, allowing the CMEs to expand more in the heliosphere, in turn allowing more CME 
detections to take place. This can also explain the reduced geoehectiveness of the CMEs. 
They have diluted magnetic energy content bec ause of their gre a ter ex pansion, and they 
interact with weaker ambient heliospheric helds: ISmith fc BaloghI (120081) earlier found that 
the radial interplanetary magnetic held component measured at high latitudes by Ulysses 
was about a third weaker during the spacecraft’s third orbit (during cycle 23/24 minimum) 
than compared to its hrst orbit (during cycle 22/23 minimum). They also found agree¬ 
ment between the Ulysees results and NASA OMNI solar wind data taken at 1 AU; the 
radial interplanetary held compone nt multiplied by the square of the rad ial distance from 
the Sun is independent of position (IBalogh et al.lll995l: ISmith et al.ll200ll) . so that in situ 
Ulysses and near-earth measurements can be compared over long periods of time. The dy- 
namical and thermal p ressures were also signihcantly smaller (by 22% and 25%, respectively 


McComas et al.l 120081) . 


This paper will compare the CME rate statistics from the three databases, SEEDS, 
CACTus and CDAW, and search for evidence of real CME rate changes since cycle 23, and 
relate these changes to decreases in the solar and heliospheric held strength. In Section [2] 
the data sets will be introdnced. The data will be analyzed and the results presented in 
Section [3l and we will conclnde in Section HI 
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Data 


LASCO has completed 18 years of nearly continuous white-light imaging of the K- 
corona, covering cycles 23 and 24. LASCO initially had three coronagraphs with overlapping 
helds of view, Cl (1.1-3 Rg), C2 (2-6 Rg), and C3 (3.7-32 Rg), where Rg is the solar radius. 
The Cl camera did not survive the temporary loss of the SoHO spacecraft in 1998. CMEs 
are therefore identihed using images from the C2 and/or C3 coronagraphs. Three o nline 
databases, the Coordinated Data Analysis Wo rkslipij/I (CDAW, Gonalswamv et ah 2009 1. the 
Solar Eruptive Events Data Systeni§ (SEEDS . Olmedo et ah 2008 1. and the Computer Aided 


CME Tracking (CACTus, Robbrecht et ah 2009 1. have been recording CMEs covering the 


LASCO era (1997-present). The CDAW CME identihcations have been achieved by visual 
inspection of the LASCO C2 and C3 coronagraph images, whereas SEEDS and CACTus 
apply automated algorithms to identify CMEs without human intervention, SEEDS using 
C2 images and CACTus using images from both C2 and C3. 

The SEEDS detection algorithm is based on projecting two-dimensional images onto 



in image sequences from LASCO C2 and C3 by constructing stackplots of C2 and C3 images, 
sharpened using the Hough transform, to detect motions of bright CME structures. Until 
2010 the algorithm ran on level 0 images from the LASCO C2 and C3 instruments and since 
2010 quick-look images from these instruments have been used. The detecti on method and 


the database are described and analyzed in detail bv iRobbrecht et al.l (1200911 . 


The CDAW database is based on lists compiled by human observers who review near- 
real-time movies of C2 and C3 images. These observers need to judge what constitutes an 
individual CME in, e.g., a close sequence of eruptions and outflows. Qualitative descriptions 
are added to many of the records in the CDAW database, such as “poor”, “very poor” and 
“C2 only”. These comments help us to determine the effects of these questionable detections 
on the overall CDAW statistics. 

We therefore have identihcations from two automated (CACTus and SEEDS) and one 
manual database, with two databases (CDAW and CACTus) based on C2 and C3 images. 


^http://cdaw.gsfc.nasa.gov 
^http://spaceweather.gmu.edu/seeds/lasco.php 
^http: //side.oma.be/cactus / 
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and one (SEEDS) based on C2 images alone. Compared to databases assembled mannally 
by human operators, automatic CME detections might be more objective as the detection 
criterion is written explicitly in a program, but the lack of manual hltering also means that 
the database needs to be treated with some caution. Qualitative comparison with the CDAW 
database is therefore a useful test of the automated algorithms’ results. The automated 
databases also provide a check of the objectivity of the human observers’ detections in the 
CDAW database. 


Comparisons between CME databases have been conducted in the past. CACTus gen¬ 
erally identihes more CMEs than CDAW, and the CACTus CME rate is better correlated 
with the sunspot cycle, thoug h the CACTus CME ra te lagged behind the sunspot cycle by 
6-12 months during cycle 23 flRobbrecht et ahl 120091) . The CDAW database is affected by 
observer bias. It has been compiled by several observers since 1997. After 2004 the CDAW 
detection rate increased signihcantly because the project began to include very narrow CMEs 
previously disregarded. 


Boursier et ahl (120091 ) found good agreement between the SEEDS and Automatic Recog¬ 


nition of Transient Events and Marseille Inventory from Synoptic maps (ARTEMIS) databases, 
perhaps reflecting the restriction of these databases to C2 images. (The ARTEMIS database 
compiled CME detections using LASCO C2 synoptic maps, but now appears to be offline.) 
The rec orded SEEDS CM E widths tend to be much narrower than the CDAW and CACTus 
widths ( Byrne et al.ll2009 ') because the SEEDS algorithm refers to C2 but not C3 images. The 
CACTus CME width distribution tends to b e scale-invariant, i n con trast to the CDAW width 
distribution which has a peak around 30° (IRobbrecht et al.N2009h . CACTus has recorded 
many more narrow CMEs than CDAW, whose human observers tend to miss many narrow 
CMEs during solar maximum. When the activity level is low, small, faint CMEs are easier to 
catch than during solar maximum. On the other hand, automated databases sometimes mis s 
wide (e.g., halo) CMEs that are picked up by human observers fe.g. iGopalswamv et al.ll2010l) . 
The observational dehnition of a CME is not uniquely dehned and varies from database to 
database, and also within a single database in the case of CDAW. 


In this paper we will revisit the issue of consistency between the three databases CAC¬ 
Tus, CDAW and SEEDS, focusing on the agreement between them regarding the question 
of changes of CME detection rates between cycles 23 and 24. These three independent solar 
eruption databases will provide a useful prohle of eruption rates during the LASCO era, that 
can be related to the solar magnetic held behavior during the same period of time. 
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3. Comparison of CME eruption distributions and rates 
3.1. Velocity distributions 
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Fig. 1.— Annual histograms of CME velocity measurements from the SEEDS database. 
Fitted log-normal functions, described by Equation ([1]), are overplotted. 


Figures [I]l3] show annual velocity histograms of detected CME velocities from the SEEDS, 
CACTus and CDAW databases, respectively. The histograms from all three databases all 
have approximately log-normal distributions, whose amplitudes have a clear solar cycle de- 
pendence. The log-norm al distribution of CME speeds has been discussed in the past by, e.g.. 


Yurchvshvn et al.l (120051) . Comparing the amplitudes of the log-normal distributions in Fig¬ 


ures m year-by-year, the CACTUS detections seem to occur in consistently lower numbers 
than the SEEDS and CDAW detections. The CDAW distributions extend to higher veloci¬ 
ties than the SEEDS and CACTUS distributions, though the CACTUS velocity histograms 
tend to peak around 300 km s“^, signihcantly higher than the SEEDS peak velocities and 
marginally higher than the CDAW peak velocities. On the other hand, the velocity distri¬ 
butions from all three databases appear to have similar widths. Figure H] shows the CDAW 
histograms with “poor” and “very poor” detections removed. The vast majority of detections 
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Fig. 2.— Annual histograms of CME velocity measurements from the CACTus database. 
Fitted log-normal functions, described by Equation ([1]), are overplotted. 

survive, and Figures [3] and 0] appear similar, with slightly reduced statistics in Figure 01 

To gain a more quantitative characterization of the velocity distributions, we model the 
overall behavior of the velocity distributions year by year using the log-normal distribution 
function. 


, , C -(logD-M^) 

y iTTav 


( 1 ) 


This function describes a peaked, asymmetric distribution and is appropriate for modeling 
the velocity distributions shown in Figures [T]01 The location of the peak, referred to as the 
mode, is given by and the median and mean by and Thus the mode < the 

median < the mean, with equality only in the symmetric case cr = 0. Overplotted on each 
histogram in Figures [I]0] is the best log-normal fit to the histogram. The free parameters of 
the log-normal function in Equation ([1]) will be used later, in Subsection 13.51 to characterize 
and compare the velocity distributions. First we examine the effects of excluding classes of 
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Fig. 3.— Annual histograms of CME velocity measurements from the CDAW database. 
Fitted log-normal functions, described by Equation ([1]), are overplotted. 

relatively unreliable detections. 


3.2. Effect of excluding narrow detections 

Figure |5] shows the changes in time of the angular width (left column of plots) and 
velocity distributions (middle and right columns of plots) associated with the CME detec¬ 
tions by SEEDS, CACTus and CDAW. The middle column of Figure |5] shows the velocity 
histograms for all CME detections and the right column shows the histograms with narrow 
cases (angular width < 30°) removed. The SEEDS data, based on C2 images only, clearly 
include a larger proportion of narrow CMEs than the other two data sets, which are based 
on data from C2 and C3. This may be because CMEs tend to expand super-radially as the 
ambient held strength drops off faster than l/(radius)^, so that as they travel from the C2 
held of view to the C3 held of view they tend to have increasing angular width. The narrow 
CMEs detected by CACTus show a strong cycle-dependence, albeit with a higher rate of 
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Fig. 4.— Annual histograms of CME velocity measurements from the CDAW database. 
Fitted log-normal functions, described by Equation ([1]), are overplotted. Poor and very poor 
detections have been excluded. 


detection since 2010. The CDAW distribution shows more of a stepwise increase in narrow 
CME detections in 2007, with a broad peak during the cycle 23 minimum in 2007-2009. 
The CACTus detection increase may be related to the LASCO C2 and C3 image cadence 
increases that occurred during 2010, whereas the CDAW peak is likely due to the inclusion 
of narrower CMEs in the database by the human observers after 2004. 


The exclusion of narrow CMEs in the right column of Figure E] has the effect of reducing 
the number of slow CMEs and moving the peaks of the log-normal distributions up the 
velocity scale. This effect is not surprising because it is well known that the ang ular widths 
and velocities of CMEs are moderately well correlated flGopalswamv et al.l 120141 1. This step 
removes most of the questionable, borderline CME detections: those with angular width 
> 30° are relatively unlikely to be time-dependent plasma outflows mistakenly identihed as 
CMEs. The velocity distributions in the right column of Figure 0 are more regular than 
those in the middle column as a result. 
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3.3. Impact of LASCO image cadence increase 

The cadence of LASCO C2 and C3 images increased by a factor of approximately two 
when other instruments on the Solar and Heliospheric Observatory (SoHO) ceased to send 
data to earth regularly, and extra bandwidth became available to LASCO. To assess the 
effect of the LASCO image cadence change on the rate of CME detection, Figure [H] shows 
a 2D histogram of annual number distributions of LASCO/C2 images per SEEDS CME 
detection and the average number of C2 images per day. Both statistics are available in the 
SEEDS database. The increase by a factor of about 2 of the nominal LASCO/C2 image 
cadence from about 50 day“^ to about 100 day“^ appears clearly in Figure O (right). This 
major change was accompanied by changes in the distribution of the number of C2 images 
per SEEDS CME detection, shown in Figure [6] (left). The figure shows that the distribution 
of images per CME detection broadened significantly when the C2 image cadence increased. 

There is no obvious increase in the maximum speed of detected CMEs between cycles 

23 and 24 in Figure O According to Figure [HI the main impact of the LASCO image cadence 
increase seems to have been to broaden the distribution of C2 images per SEEDS CME 
detection rather than introducing more detections of very fast or small, faint CMEs only 
detectable in a couple of high-cadence images. If the leading effect of the C2 image cadence 
increase had been to introduce many more detections then these would have included a large 
increase in detections in the left bins of the left plot of Figure [HI either from very fast CMEs 
caught in just a few high-cadence images, or small, faint CMEs, caught in just a few images 
before becoming undetectably diffuse. However, these bins are less populated during cycle 

24 than during cycle 23, a side-effect of the broadening of the distributions. 


3.4. Effect of excluding “poor” and “very poor” CDAW detections 

Figure [7| shows 2D histograms of the CDAW CME velocities and angular widths for the 
cases including all detections, and with “poor” and “very poor” detections removed. The 
right column of plots shows the velocity distributions for cases with angular width > 30°. 
The top row of Figure [7| is identical to the bottom row of Figure [HI The curious maximum 
in the number of narrow CMEs coinciding with solar activity minimum (2007-2009) is again 
clearly visible in the left plots, including the plots with “poor” and “very poor” detections 
removed. To avoid artifacts associated with observer bias it is evidently not sufficient to 
remove the “poor” and “very poor” detections. Implausibly large post-2004 statistics are 
evident in the columns up to 30° in all three left plots of Figure [3 Beyond 30° the histograms 
resemble the CACTus histogram in Figure [H] (middle left). The main difference between the 
left plots of Figure [3 for CMEs with angular width > 30° is that the post-2004 statistics are 
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clearly larger in the top left plot than the pre-2004 statistics because of the observer bias. 
This artifact is effectively removed by the exclusion of “poor” and “very poor” detections. 
We can verify this by qualitatively comparing the bottom left plot of Figure [7] to the middle 
left plot of Figure E] which shows the distribution from the more objective CACTus algorithm 
processing the same C2 and C3 images. 

The velocity distributions in the middle column of Figure [7] include detections of all 
angular widths. For comparison, the equivalent velocity distributions for detections with 
angular width > 30° are plotted in the right column of Figure [71 As in Figure El the main 
effect of excluding narrow CMEs is to reduce the number of slow CMEs. Comparing the 
plots in the middle and right columns of Figure [71 removal of the “poor” and “very poor” 
detections seems to have only a minor effect on the velocity distribution, whereas removing 
narrow CMEs has a major simplifying effect, eliminating irregularities in the distribution 
and bringing it into much closer qualitative agreement with the CACTus distribution in 
Figure El 


Figures [71 and El indicate that we can make an improved qualitative comparison between 
CACTus and CDAW statistics if we conhne our attention to CMEs with angular width 
> 30° and ignore “poor” and “very poor” CDAW detections. Figure [71 shows also that the 
statistics of halo CME detections (angular width 360°) are not changed much by the exclusion 
of “poor” and “very po or” detections. Halo CME s form a small but notably robust subset 
of the CDAW database flCopalswamv et ahl [20151) . 


3.5. Comparison of the log-normal distribution parameters in time 

It is apparent from Figures El and [71 that the maximum speed of detected CMEs did not 
signihcantly increase as a result of the cadence change, nor is there an obvious increase in 
the median velocity. We can address the question of the median and mean velocities more 
qualitatively using Figure [HI where the prohles of the modes, medians and means of the log¬ 
normal best hts to the velocity distributions are plotted. These best-£t log-normal functions 
are derived by optimizing the free parameters of Equation Narrow CMEs (angular 
width < 30°) are excluded from these calculations. The resulting time-prohles for the modes, 
medians and means of the log-normal best-hts shown in Figure El are strikingly similar for the 
CDAW and CACTus data, but somewhat different for the SEEDS data. The resemblance 
of the CDAW and CACTus plots to each other is likely due to the fact that the CDAW and 
CACTus databases are both based on C2 and C3 data, and the removal of the narrow CME 
detections has eliminated most of the effects of subjective CME identification in the CDAW 
database. The CDAW plot shows three sets of parameter prohles for CME detections with 
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angular width > 30°: all such detections including “poor” and “very poor” detections (black 
lines), with “poor” detections excluded (dark grey lines), and with “very poor” detections 
excluded (light grey lines). The exclusion of “poor” and “very poor” detections does not 
produce a dramatic effect on the log-normal velocity distributions, but it increases the mean 
velocity by up to 10% during the late maximum and decline of cycle 23 and cycle 24. This 
change is much smaller than the shift from faster average velocities of about 550 km s“^ at 
the peak of cycle 23 to a value of about 400 km s“^ typical for cycle 24 so far. A similar 
difference between cycle 23 and 24 peaks is seen in the CACTus data supporting the notion 
that a real physical change in velocity occurred between the two cycles. However, it appears 
that the enhanced cycle 23 peak only occurred in the two databases that extended to the 
C3 held of view: CACTUS and CDAW. The SEEDS mean and median velocities have only 
modest peaks in 2005, as shown in Figure [HI Otherwise the SEEDS parameters maintain 
steady values throughout the active years of cycles 23 and 24. In particular, there is no 
obvious change in the velocity distributions between late cycle 23 (2004, 2005) and early 
cycle 24 (2011-2014). 

The mode, median and mean values associated with the log-normal hts to the SEEDS 
statistics for C2 images CME“^ (FigureEl left) are shown in the bottom right plot of FigurelHl 
The mean and median C2 images CME“^ show a clear stepwise change coinciding with the 
C2 image cadence change in 2010. The SEEDS velocity parameters do not show such a 
stepwise change. The CDAW and CACTus velocity parameters also have similar values 
before (2004, 2005) and after (2011-2014) the cycle 23 minimum. The agreement of all three 
databases on this point suggests that the LASCO image cadence change did not signihcantly 
impact the velocity distributions of the CMEs represented in Figure [HI those with angular 
width > 30°. 


3.6. Comparison of CME rates, and the effects of excluding narrow and 

poorly detected CMEs 

Figure [9] shows the annual CME rates for the SEEDS, CACTus and CDAW databases, 
including all detections (left plots) and for detections with angular width > 30° (middle 
plots). Also shown are the amplitudes of the log-normal function hts for detections with an¬ 
gular width > 30° (right plots). The agreement between the middle and right plots indicates 
that the log-normal functions represent the observed velocity distributions satisfactorily. In 
each plot the annual average monthly sunspot number is over-plotted, scaled according to 
the linear regression coefficients, and standard deviations of annual means are indicated both 
for CME rates and the sunspot number. This allows us to track the deviation of the CME 
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rates from the sunspot number. Whereas in the p ast the CME rate has been well corre¬ 
lated with the sunspot number (e.g., during cycle 2l lWebb fc Howardlll994l) . it is clear from 
Figure in] that these LASCO-based CME rates deviate from the sunspot numbers in a sta¬ 
tistically signihcant manner. In all cases the CME rate is systematically higher per sunspot 
number during cycle 24 than during cycle 23. On the other hand, the CME rates from the 
different databases differ signihcantly from each other in general. The CME detection rates 
evidently depend on differences in how the databases are constructed, including whether 
manual or automated detection methods were used, and whether the detections were based 
on C2 and/or C3 images. 


As we saw in Figure [T] Figure [9] also shows that the manually-derived CDAW statistics 
are affected by the inclusion of narrow CMEs in 2005 and some observer-dependent artifacts, 
giving the cycle 23 minimum disproportionately numerous detections and making the cycle 
24 maximum peak much higher than the cycle 23 maximum. Again we see here, comparing 
the CDAW plots to the CACTus plots, that these artifacts are largely removed by excluding 
the CMEs with angular width < 30°. 


Also shown in Figure IH] are CDAW CME rates excluding “poor” and “very poor” de¬ 
tections. The effects of these exclusions are more modest, but they do dampen the cycle 24 
maximum peak so that it is approximately the same height as the cycle 23 maximum. The 
cycle 23 and 24 maximum peaks are also approximately the same height in the CACTus 
plots, including or excluding the CMEs with angular width < 30°. 


In the plots of the SEEDS CME rates in Figure [H the cycle 23 and 24 maxima are 
approximately the same height when all detections are included, but the exclusion of narrow 
(angular width < 30°) CMEs shrinks both peaks by more than a half, and the cycle 24 
peak becomes about 10-20% smaller than the cycle 23 peak. Nevertheless there is a surplus 
of CMEs relative to the sunspot number after the cycle 23 sunspot maximum (2000-2002) 
compared to before. This surplus is reduced by the exclusion of narrow CMEs, those with 
with angular width < 30°. However, the surplus of CMEs per sunspot number in cycle 24 
compared to cycle 23 remains large and signihcant for CMEs with for angular width > 30°. 
This signihcant enhancement of CME detections per sunspot number for cycle 24 relative to 
cycle 23 is therefore common to the SEEDS, CACTus and CDAW data. 


The CACTus and CDAW curves in Figure [9]for angular width > 30° also show signihcant 
increases relative to the sunspot number during the decline of cycle 23, in 2004 and 2005. 
Since this surplus appears in the CACTus and CDAW data, it seems to be independent of the 
observer biases in the CDAW data. However, thi s increase is absent from the SEEDS plot in 
Figure [9] for angular widths > 30°. According to lLamv et ahl fj201411 . ARTEMIS recorded an 
enhanced CME rate during the rise of cycle 24 but no signihcant enhancement during the 
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decline of cycle 23. This would agree with the SEEDS > 30° plot in Figure[9l Such agreement 
may be related to the exclusive use of C2 images by both SEEDS and ARTEMIS. That this 
late cycle 23 surplus occurs in the CACTus and CDAW data but not so prominently in the 
SEEDS (or ARTEMIS) data suggests that these extra CMEs may have required C3 images 
for their detection. 

We now look at the effect of the LASCO image cadence change on the marginal CDAW 
detections, and the impact of these marginal detections on the overall statistics. Figure ITU] 
shows the average daily rate of detections by CDAW where the CME appeared in only 2 
or 3 LASCO images, or only in C2 images, or only in C3 images. These rates are plotted 
separately for the cases with all angular widths and excluding those whose angular widths 
exceed 30°. Poor and very poor events are excluded. All of the rates plotted in Figure dU] 
are low, and are much reduced in cycle 24 compared to cycle 23. The rate of CMEs detected 
in only 2 or 3 images was about one per 10 days (one per 20 days for angular widths > 30°) 
during cycle 23, and fell to zero in cycle 24. C2-only detections occurred at a rate of about 
one every 3-4 days during cycle 23 (one every 5-6 days for angular widths > 30°) and fell 
to one per 10 days in cycle 24. C3-only detections occurred every 15 days or so on average 
during cycle 23, almost all with angular width > 30°, and fell to one every 30 days in cycle 
24. The LASCO image cadence changes for cycle 24 therefore seem to have signihcantly 
affected these marginal CME detections. As we discussed in the context of Figure O the 
main effect of LASCO image cadence change seems to have been to increase the number 
of images per CME detection and to decrease the number of marginal detections as shown 
in Figure [TUI However, their impact on the total CME statistics in Figure 0 was evidently 
small. 


3.7. The cycle 23/24 CME rate increase, and the decrease in the solar and 

heliospheric magnetic field strength 

The top panel of Figure [TT] shows a combined plot of the CME rates from SEEDS, 
CACTus and CDAW, for detections with angular width > 30°. Poor and very poor detections 
are excluded in the case of CDAW. The sunspot number is also over-plotted, and the CME 
rates and sunspot number are all normalized by their average values during the cycle 23 
maximum years 2000-2002. Thus, the divergence of the CME rates from the sunspot number 
since the cycle 23 maximum is emphasized. In 2004, during the decline of cycle 23, the 
CACTus and CDAW CME rates diverged sharply from the sunspot number. Ever since 2004 
these CME rates have remained elevated relative to the sunspot number by a statistically 
signihcant margin. Curiously, the SEEDS CME rate did not diverge from the sunspot 
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number until the onset of cycle 23 in 2010-2011. The SEEDS CME rate caught up with the 
CACTus and CDAW rates in 2011, and has remained signihcantly elevated relative to the 
sunspot number throughout cycle 24 so far. 

The second panel of Figure ITT] plots the mean photospheric flux density from 1997 to 
2015. The photospheric held is generally well correlated with the sunspot number, being 
dominated by active region hux during activity maxima. The photospheric hux has been 
around 40% lower since the decline of cycle 23 compared to the ascent of cycle 23. The polar 
and equatorial dipole components of the photospheric held are plotted in the third panel of 
Figure [TTJ The polar dipole represents the strength of the polar helds while the equatorial 
dipole follows the activity cycle. Both dipole components rehect the 40% decrease in the 
photospheric held that occurred during the decline of cycle 23, around 2003-2004. 


The fourth panel of Figure dH plots the OMNI2 radial IMF cor nponent at 1 AU. The 
radial IMF has been widely reported to have decreased by about 30% fjSmith fc Baloghll2008[) 
between the cycle 22/23 and 23/24 minima, but the plot shows that this decrease actually 
occurred during the decline of cycle 23, early in 2004. During the cycle 23 declining phase, 
the photospheric polar helds ceased to strengthen around 20 03, even thoug h there was much 
magnetic activity between 2003 and the end of the cycle. iPetriel (120121 ) argued that this 
was due to the the positive and negative active region hux latitude centroids statistically 
coinciding in each hemisphere, i.e., the Joy’s law dipole tilt becoming insignihcant, from 
around 2003 until the end of the cycle. This left the polar helds about 40% weaker after 
their cycle 23 polarity reversal compared to before. This weakening of the photospheric held 
during the decline of cycle 23 produced a nearly stepwise 30% decrease in the radial IMF in 
early 2004. On the rise of cycle 24 in 2010, the IMF only rose to cycle 23 declining-phase 
levels, about 30% short of cycle 23 maximum-phase values. Recently the IMF has spiked as 
the polar helds have strengthened before the activity has decayed away. 


It would be useful to relate the photospheric held to the IMF via a magnetic model for 
the coronal op en hux. The obvious way to achieve this is the p otential-held source-surface 
(PFSS) model f Schatten et ah 1969 : Altschuler fc Newkirk 1969 ). The solar open hux should 
match the radial component of the IMF integrated over the sphere at 1 AU. However, the 
PFSS model cannot be relied upon to model the open hux over full solar cycles with a hxed 
outer bou n dary (source surface) radius. Using the source surface radius as a free parameter. 


Lee et ahl (120091) showed that PFSS coronal hole maps and open huxes better reproduced 


observed coronal hole distributions and IMF measurements with smaller source surface radii 
during the cycle 24 minimum than during the cycle 23 minimum, sugges ting that the s ource- 
surface radius should be lowered when the photospheric held is weak. lArden et al.l (120141) 
found that the source-surface radius “breathes”, such that its height needs to be around 
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2.5 solar radii during solar minima and 15-30% higher during solar maxima, according to 
comparisons between PFSS open fluxes and the measured IMF. This situation is further com¬ 
plicated by the general lack of quantitativ e agreement betw een photospheric magnetic held 
measurements from different observatories flRilev et al.ll2014l) . We defer further investigation 
of the link between the photospheric and interplanetary fields to a future paper. 

Comparing the panels of Figure fTTl the timing of the CACTus and CDAW CME rates’ 
divergence from the sunspot number, in 2004, matches the decrease in the IMF. A causal 
physical link bet ween the decreased I MF a nd an increased CME rate per sunspot number 
was suggested bv iGopalswamv et al.l fj20141 ). who argued that a decreased heliospheric total 
(magnetic-fplasma) pressure has allowed CMEs to expand more as they travel from the solar 
surface into the heliosphere, thereby enhancing the CME rate statistics. 

The energy density of the radial magnetic held (that produces the transverse magnetic 
pressure stresses that woul d act upon an expanding CME) fall s off much less rapidly than 
that of the transverse held (jSchatten et al.lll969l: ISchattenlll97lh . The energy density of the 
transverse held approximately balances the plasma pressure at about 0.6 solar radii above 
the photosphere. Thus the plasma extends the magnetic held outward near this point. In 
the case of the radial held, equal ity with the pl asma energy density is only reached at the 
Alfven point near 25 solar radii flSchattenlll97l[) . The magnetic held is therefore expected 
to constrain the plasma expansion out to to 25 solar radii, covering almost all of the fields 
of view of the LASCO C2 and C3 coronagraphs. Although IGopalswamv et al.l (120141) were 
comparing CME statistics from the rise phases of cycles 23 and 24, their argument links the 
CME rate to the solar and interplanetary magnetic held whose decrease in strength occurred 
during the decline of cycle 23, around 2003-2004. If the argument is correct then one would 
expect the increase in CME rate per sunspot number to begin around 2004, as the CACTus 
and CDAW statistics clearly do in the top panel of Figure HH 

The question remains why the CACTus and CDAW CME rates diverged from the 
sunspot number in 2004, whereas the SEEDS rate waited until 2010-2011 before diverg¬ 
ing. It is difficult to explain this in terms of a changing rate of CME expansion between 
the C2 and C3 fields of view. Such an explanation may again require improved coronal 
field modeling to estimate the magnetic (and total) pressure change between the C2 and 
C3 fields of view over time. It is worth noting, however, that the Nobeyama radio helio¬ 
graph prominence eruption statistics also show a statistically significant in eruption rate per 
sunspot number from cycle 23 to cycle 24 , and also a smaller but statistically significant 
increase per sunspot number around 2 004 (^triel 201^. Because eruptive prominences are 
almost always accompanied by CMEs (iMunro et al.lll979r). prominen ce eruptions and CMEs 
can often be identified with each other (IGopalswamv et al.l l2003bl) . hence this increase in 
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prominence eruptions around 2004 supports the notion that the CME rate also rose at these 
times. 


4. Conclusion 


The results are summarized as follows. In the three LASCO CME databases SEEDS, 
CACTus and CDAW, a statistically signihcant increase in the rate of CME detections with 
angular width > 30° per sunspot number was found for cycle 24 compared to cycle 23. In the 
two databases based on both LASCO/C2 and C3 images, CACTus and CDAW, the upward 
divergence of the CME rate relative to the sunspot number began in 2004 after the polar field 
reversal. At nearly t he same time, the IMF decreased by ~ 30%. These results are consistent 
with the the result of iGopalswamv et al.l (12014 l2015l ) linking enhanced halo CME detections 
to increased CME expansion in a heliosphere of decreased total (magnetic+plasma) pressure. 
On the other hand, the SEEDS CME rate did not diverge from the sunspot number until 
the rise of cycle 24, in 2010-2011. It is possible that this can be explained by the restriction 
of the SEEDS detections to the C2 held of view. Such an explanation would need to rely on 
an improved model for the cycle dependence of the global coronal held strength at 2-3 solar 
radii than is currently available. 


The LASCO C2 and C3 image cadence changes in 2010 may have some ehect on the 
CME detection rates. However, this ehect is likely to have been small in view of the increase 
in the number of images per detection, the decrease in low-quality detections and the lack of 
evidence of enhanced detections of very fast or faint CMEs only detectable in high-cadence 
sequences of images. Our restriction to angular widths > 30° focuses our study almost 
exclusively on CMEs that would be equally well detected using cadences of 10 or 20 minutes. 
The increases in CME detection rate increases seem too large to be explained by marginal 
detection changes, and the explanation in terms of enhanced CME expansion in a coronal 
medium of reduced transverse magnetic pressure beginning in 2004 is consistent with the 
CACTus and CDAW data. The later increase of the SEEDS CME rate, at the beginning of 
cycle 24, remains to be explained. 


The author thanks the referee for helpful comments. SOHO is a project of international 
cooperation between ESA and NASA. The CDAW CME catalog is generated and maintained 
at the CDAW Data Center by NASA and The Catholic University of America in cooperation 
with the Naval Research Laboratory.This paper uses data from the CACTus CME database, 
generated and maintained by the SIDC at the Royal Observatory of Belgium. The paper also 
uses CME statistics from the SEEDS project at George Mason University’s Space Weather 
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Laboratory, that has been supported by NASA’s Living With a Star Program and NASA’s 
Applied Information Systems Research Program. SOLIS data used here are produced co¬ 
operatively by NSF/NSO and NASA/LWS. NSO/Kitt Peak SPMG data used here were 
produced cooperatively by NSF/NOAO, NASA/GSFC, and NOAA/SEL. The OMNI data 
were obtained from the GSFG/SPDF OMNIWeb interface at http://omniweb.gsfc.nasa.gov 
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Fig. 5.— Stacked annual histograms of CME angular widths (left) and angular widths for 
all detections (middle) and for detections with angular width > 30° (right), from the SEEDS 
(top), CACTus (middle) and CDAW (bottom) databases. The dotted lines in the left plots 
indicate the angular width cutoff at 30° below which the detections are excluded from the 
right plots. 
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Fig. 6.— Stacked annual histograms of the number of LASCO C2 images per detected CME 
in the SEEDS database (left) and of the number of C2 images per day (right). 
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Fig. 7.— Stacked annual histograms of CME angular widths (left) and angular widths for 
all detections (middle) and for detections with angular width > 30° (right), from the CDAW 
database. All detections are included in the top plots, very poor detections are excluded 
from the middle plots, and poor and very poor detections are excluded from the bottom 
plots. The dotted lines in the left plots indicate the angular width cutoff at 30° below which 
the detections are excluded from the right plots. 
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Fig. 8.— The annual mean (solid lines), median (dashed lines) and mode (dotted lines) of 
the log-normal velocity distributions in Figures [T]ll] are plotted against time for SEEDS (top 
left), CACTUS (top right) and CDAW (bottom left) detections. Here cases with angular 
width < 30° have been excluded from all databases. Eor CDAW data, the different shades 
of the curves represent calculations including “poor” and “very poor” detections (black), 
excluding “very poor” detections (dark grey) and excluding both “poor” and “very poor” 
detections (light grey). The bottom right plot shows the annual mean (solid lines), median 
(dashed lines) and mode (dotted lines) of the log-normal distributions of the annual SEEDS 
statistics for the number of LASCO C2 images per detected CME as functions of time. 
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Fig. 9.— Average daily CME rates from SEEDS (top row), CACTus (second row) and 
CDAW (3rd-5th rows) with all cases included (3rd row), very poor cases excluded (4th row) 
and poor and very poor cases excluded (5th row). In the hrst column, CMEs of all angular 
widths are included, whereas in the middle column only CMEs with angular width > 30° are 
included. The error bars indicate the standard deviations of the annual means. The right 
column shows the amplitude of Equation ([1]) htted to the annual distribution of detections 
with angular width > 30°. The average monthly sunspot number is overplotted in dashed 
curves for comparison. 
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Fig. 10.— Average daily CME rates from CDAW where the CME appeared in only 2 or 3 
LASCO images (left), only in C2 images (middle), and only in C3 images (right). Poor and 
very poor events are not included. Solid lines show the rates for all CMEs and dashed lines 
rates for CMEs with angular width > 30°. 
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Fig. 11.— Normalized annual average daily CME rates (top panel) from SEEDS, CACTus 
and CDAW including only detections with angular width > 30°, for the years 1997-2014. All 
poor and very poor cases have been excluded from the CDAW data. The normalized average 
monthly sunspot number is over-plotted in dashed curves for comparison. The CME rates 
have been normalized by their average values during the cycle 23 maximum years 2000-2002. 
Thus, the plot shows the divergence of their evolution since the cycle 23 maximum. The 
second panel shows the total photospheric magnetic flux from NSO KPVT and SOLIS/VSM 
synoptic magnetograms. The third panel shows the polar and equatorial dipole components 
of the same synoptic magneto grams. The bottom panel shows the OMNI 2 radial magnetic 
field component at 1 AU. 








